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The amplitude ratio of the susceptibihty (or second size-moment) for two-dimensional percolation 
is calculated by two series methods and also by Monte-Carlo simulation. The first series method 
is a new approach based upon integrating approximations to the scaling function. The second 
series method directly uses low- and high-density series expansions of the susceptibility, going to 
unprecedented orders for both bond and site percolation on the square lattice. Putting all methods 
together we find a consistent value r~/r"'" = 162.5 ± 2, a significant improvement over previous 
results that placed the value of this ratio variously in the range of 14 to 220. 



Percolation is one of the fundamental problems in sta- 
tistical mechanics 0, 01 and is perhaps the simplest sys- 
tem exhibiting critical behavior. Through its mapping 
to the g-state Potts model (for q ^ 1) many theoret- 
ical predictions follow, such as exact critical exponents 
in two dimensions. Yet many unanswered questions re- 
main. One of these is the value of amplitude ratios, which 
represent universal ratios of quantities related to inte- 
grals of the scaling function. A great amount of work 
has been done investigating amplitude ratios of various 
systems, both to demonstrate that systems expected to 
be in a given universality class have the same ratios, and 
to determine their values accurately 0. The study of 
amplitude ratios remains an active area of research (e.g., 

iiiHiiMini). 

In this paper we study specifically the universal ampli- 
tude ratio r^/r+ for percolation, where T~ and r+ are 
the amplitudes of the second size-moment (also called 
the susceptibility) in the low- and high-density phases, 
respectively. This ratio has been especially difficult to es- 
timate accurately and the status of the estimates is very 
controversial. In |^ values are quoted in a wide range 
from 14 to 220 based on numerical estimates from Monte- 
Carlo simulations and series expansions for various per- 
colation models. Furthermore, there are no reliable field- 
theoretical estimate for this quantity. Several years ago 
Dclfino and Cardy [l3| studied the g-state Potts model 
using methods from quantum field theory and predicted a 
value of 74.2 for percolation, by extrapolating the results 
for q = 2, 3, and 4 to g = 1. This was consistent with 
the numerical work of Corsten et al. [l^, who gave the 
value 75 (-1-40,-26), but inconsistent with many other 
measurements In this work, we study bond and site 
percolation on the square lattice and use extensive exact 
enumerations to obtain estimates for r~/r"'", using two 
different approaches: one a novel approach based upon 
directly integrating approximations to the scaling func- 
tion, and the second a more conventional analysis of the 
high- and low-density series for the susceptibility. The 
estimates for the amplitude ratio are consistent with the 
value r~/r+ — 162 ± 3. We also carried out a Monte- 



Carlo calculation, which gave an almost identical value 
of 163 ± 2. These results are a significant improvement 
on previous published numerical estimates. 

Percolation models are commonly formulated in a lat- 
tice setting with the edges and/or vertices occupied (or 
vacant) with probability p { or 1 — p). In this paper we 
limit our study to bond and site percolation on the square 
lattice Z^. We shall refer to occupied edges and vertices 
as bonds and sites, respectively. Nearest neighbor bonds 
(sites) are said to be connected and clusters are sets of 
connected bonds (sites). The behavior of the model is 
controlled by the occupation probability (or density of 
bonds/sites) p. When p is smaller than a critical value 
Pc all clusters remain finite. Above Pc there is a non-zero 
probability of finding an infinite cluster. The critical oc- 
cupation probability is know exactly for bond percola- 
tion, Pc = 1/2 14], and to a high degree of numerical 
accuracy for site percolation, pc ~ 0.59274621(13) flU . 
The average cluster size S{p) T^{pc — p)~'^ , which di- 
verges as p — > , plays a role similar to a susceptibility, 
with 7 = 43/18. 

Percolation problems are closely related to the combi- 
natorial problem of the enumeration of lattice animals, 
which are connected subgraphs of a lattice. The size 
of a lattice animal is the number of connected sites (or 
bonds). A vertex (or edge) is said to be a perimeter site 
if it is a nearest neighbor of a site in the lattice animal. 
Series expansions for various percolation properties, such 
as the average cluster size, can be obtained as weighted 
sums over the number of lattice animals, gs,ti enumerated 
according to the number of sites (bonds) s and perimeter 
t. Perimeter polynomials are defined as 

DM^Y.9s.tq'- (1) 

t 

We have calculated the perimeter polynomials up to size 
35 (bond) and 40 (site). The central property describing 
the cluster statistics in percolation is rij, defined as the 
number of clusters (per site) containing s occupied sites 
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or bonds, as a function of the occupation probabihty p, 

nM^P'Ds{l-p) = J2gs,tP%l-pY. (2) 
t 

The scahng hypothesis states that Ug behaves as 

ns{p) ^ cqs^^ f{ci{p ~ pc)s'') {p ^ Pc, s oo), (3) 

where the critical exponents r and a and the scaling func- 
tion / are universal, while co and Ci are non-universal 
metric factors. In two dimensions t — 187/91 and 
a = 36/91. We shall often use the scaling variable 
z — ci{p ~ Pc)s" ■ It then follows that 

S{p) = ^s^ns{p) J cos'^^''f{ci{p-Pc)s'')ds 



a 
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\p-Pc\-'' \ izr^+v(^)dz 



(4) 



where 7 = (3 — t)/(j. The integration in the high- or 
low-density case extends from to 00 or from —00 to 
0, respectively. The amplitudes F"*" and T~ are given 
by the non-universal constant cq(\^ j a multiplied by the 
corresponding integrals of the universal scaling function 
/(z). It thus follows that the ratio r~/r+ is universal. 
Now, from the knowledge of ns{p) for finite s, where Usij)) 
is a polynomial in p, we can estimate the integrals in (0J) 
by approximating f{z) by ns(z) = s'^ns(z), where z — 
{p — Pc)s'^ , and taking the ratio of the relevant integrals 
involving ns{z). Obviously, ns{z) is just a polynomial 
approximation to f{z) so we can't extend the integration 
to infinity (the integral would just diverge). However, 
there is a natural cut-off provided by the scaling variable 
z and the fact that the physical low-density region is 
< p < Pc, so the integral involving z over the low- 
density region runs over the interval [—z^ — —s'^pcO]. 
Likewise, integrals over the high-density region include 
the interval [0, s°'(l — Pc) = The integrals are easily 
evaluated since ns{z) — OfcZ*^ are polynomials, which 
we can determine by exact enumeration. We get in the 
high-density case: 



1+ ^ / %p-in,(z)dz = f z''-^^akz''dz 

Jo Jo fc^o 

= ^afe(z+)''^+V(fc + 7), 



k=0 



and in the low-density region 

= [\~zy~^ns{z)dz^ rz''-i^afc(-z)'=d2 
J-z^ Jo fc^o 

^ ^(-l)'=afe(z_)'=+V(fc + 7), 



from which we obtain the estimate rj/r+ = I~ /I^. We 
shall call this ratio rs for short. The amplitude ratio 
r^/r+ is obtained from the limit of as s —^ 00. 

In Fig.^we show (in the left panel) the estimates of the 
ratio rs — I'^ /It for bond percolation. The estimates 
display some curvature when plotted against 1/s, but 
seem to become close to a straight line and extrapolate 
to a value around 160. In order to extrapolate to the 
limit s — > 00 with a bit more confidence we need to try 
and work out the asymptotic form of rg- In the middle 
panel we show a log-log plot of the difference between 
consecutive ratios, dg — rg — rg-i, against 1/s. Clearly, 
ds has a power-law decay with 1/s. This means that 
rs « r^/F"*" -I- a/s", and is is straightforward to show 
that then ds oc l/s"+^. We now try to estimate the 
decay exponent a. In the right panel of Fig. ^ we plot 
estimates for the local exponent as against 1/s, where 
as is obtained from a linear regression of logdj vs. logs 
using the 5 values from s to s — 4. From this figure we 
estimate that a -I- 1 = 1.85(5). 

Next we extrapolate the data for by fitting directly 
to an assumed asymptotic form 



r+ 



(5) 



where the exponents at form a strictly increasing se- 
quence. By way of justification we can mention that 
similar forms are commonly found in the study of the 
asymptotic behavior of series coefhcients and arise from 
corrections to scaling. We don't know the values of a^, 
except for the leading exponent a^ = a, so we assume 
that ai ^ a + i (this is akin to including only one non- 
analytic correction to scaling), and fit to the form 



F- 



(6) 



we take a sub-sequence of terms 
Ts-k}, plug into the formula above 



fe=o 



That is 

{rs,rs~i, 

and solve the fc -I- 1 linear equations to obtain estimates 
for the amplitudes. It is then advantageous to plot 
estimates for the leading amplitude F~/F+ against 1/s 
for several values of fc. The results using the value 
a = 0.85 are plotted in Fig. |21 We notice that the 
estimates obtained with fc = 2, 3, and 4 are quite stable 
and show little dependence on s for large values. This 
would indicate that rs is well approximated by the 
assumed asymptotic form. We estimate from this that 
F~/F+ = 159.2 ± 0.2, where the error bars represent 
fluctuations but not systematic error associated with 
this method. We also tested the sensitivity of the ex- 
trapolation procedure to the value of a by fitting to the 
form © using the values a = 0.8 and 0.9, respectively. 
The resulting estimates for the amplitude ratio were 
only marginally lower or higher and the procedure is 
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FIG. 1: The amplitude ratio for bond percolation (left panel), a log-log plot of the difference between consecutive ratios 
(middle panel), and the local exponent of the difference plot (right panel). 
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FIG. 2: Estimates of the amplitude ratio P /F"*" for bond 
percolation from fits to the form (|SJ with a — 0.85. 



not very sensitive to changes in a, at least within a 
reasonable range. 

A similar analysis was carried out for site percolation. 
In this case the estimates for rg displayed a pronounced 
curvature when plotted against 1/s and the extrapolation 
to the limit s — > oo was more challenging than in the 
bond case. Wc estimate that r^/r+ = 164.5 ± 1.5. 

For our second approach, we have obtained estimates 
for the amplitudes T~ and F"*" directly from the low- and 
high-density series for the second size- moment M2{p). 
In the low-density region the series was simply obtained 
from the perimeter polynomials 

M2(p) = J2 E(i - - r-(p, -py, 

s t 

while the high-density series was calculated separately to 
order 51 (bond) and 55 (site). Estimates for the ampli- 
tudes were obtained by transforming the original series, 
using our knowledge of the critical point and exponents, 
into series which have a simple pole at Pc, with a residue 
from which we can calculate the amplitude. The fol- 
lowing two methods are quite standard [l6l |. We can 



raise the series to the power I/7 to get a series with 
the behavior {T~Y^'' /{pc — p) or we can look at the se- 
ries M2{p){pc — p)"'~^ which should have the behavior 
T~/{pc —p)- The amplitudes can then be estimated by 
forming Fade approximants to the transformed series and 
calculating the residues. Another approach is to get com- 
pletely rid of the physical singularity at pc by studying 
the transformed series M2(j)){pc—pV oi" M^2{pY^^ {pc—p) 
and evaluating Fade approximants to this series at pc- 
As pointed out by Daboul et al. in their recent study of 
2D percolation [lj|, such biased methods for calculating 
the amplitudes generally lead to more accurate estimates 
when judged purely on the spread among the estimates 
as obtained from individual approximants. However, this 
higher apparent accuracy can be misleading in that bi- 
ased approximations can display systematic corrections 
on a scale exceeding the fluctuations. Thus great care 
must be taken and one should be careful not to rely 
too heavily on the biased estimates, particularly when 
it comes to determining the error bars on the estimates. 

In Fig. 13 we have plotted estimates for the high- 
density amplitude F for site percolation. The esti- 
mates are obtained from Fade approximants to the series 
jsM2{q){qc - q)^~'^ (upper panel) and jrM2{q){qc - qV 
(lower panel), where the factor 1/q'* is included in or- 
der that the transformed series has a non-zero constant 
term. As expected the biased estimates (lower panel) 
are quite well-converged and appear to settle down at 
a value around 0.02208(2). This is obviously slightly at 
odds with the unbiased estimates (upper panel) which 
favor a higher value, but these estimates have a greater 
spread and a pronounced downwards drift. Estimates 
from the other methods display similar trends and in 
particular have a greater spread than those in the lower 
panel. To be cautious we adopt the rather conserva- 
tive estimate F = 0.02215(15) and from this we get 
F+ = q^f = 0.000609(4). Similarly we estimate that in 
the low-density case F" = 0.09819(6) and thus r-/r+ = 
161.2 ± 1.2. From a similar analysis we estimate that 
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at the 68% confidence interval. We have thus found con- 
sistent precise series-expansion and Monte-Carlo results, 
laying to rest the controversy on the value of this am- 
plitude ratio. Most likely, the wide spread in previous 
measurements were due to the large finite-size effects in 
Monte-Carlo simulations, and the relatively short series 
used in exact enumerations. We have eliminated both of 
those shortcomings in the present work. 

As part of a larger on-going project we have also cal- 
culated perimeter polynomials on the honeycomb lat- 
tice, and a preliminary analysis yields the estimates 
r-/r+ = 166 ± 5 (bond) and 170 ± 7 (site). Finally, 
we mention in passing that we have studied the area mo- 
ments of bond percolation on the square lattice (where 
the area is the volume enclosed by the external hull walk) 
(see [l^ for details) and find that the amplitude ratio of 
the second area-moment is 175 ± 10. This value is close 
to the one for the second size-moments, but it is expected 
to have a somewhat different value because it represents 
a different moment of the scaling function f{z). 

This research has been supported by the Australian 
Research Council (IJ) and the NSF (RMZ) under the 
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paper used the computational resources of the Australian 
Partnership for Advanced Computing (APAC) and the 
Victorian Partnership for Advanced Computing (VPAC). 



FIG. 3; The high-density amplitude F for site percolation. 



the amplitudes of the second size-moments of bond per- 
colation in the low- and high-density phases are F^ = 
0.15001(8) and F+ = 0.0009290(15), which yields the es- 
timate F^/F+ = 161.5 ± 0.4 for the amplitude ratio. 

The estimates from the series studies are consistent 
with the value r~/r+ = 162 ± 3, which is a significant 
improvement over previous results. We have also carried 
out a Monte-Carlo test of this amplitude ratio, study- 
ing bond percolation on the square lattice. The method 
used was to generate individual clusters by a growth al- 
gorithm on a large lattice at p = 0.47, 0.48, 0.49, 0.495, 
0.505, 0.51, 0.53, and 0.53. For p < p^, all the clus- 
ters terminate, while for p > pc, the clusters that keep 
growing beyond a clear cutoff can be identified as being 
part of the infinite cluster and discarded. This method 
proved superior to previous MC works which generally 
count clusters on fully populated lattices and have sig- 
nificant finite-size effects; here, there were essentially no 
finite-size effects as long as p is kept sufficiently away 
from Pc and the lattice is made large enough. This work 
yielded the result F-/F+ = 163 ± 2, which was quoted 
previously in [l8| . 

The two determinations of the amplitude ratio (by se- 
ries and Monte-Carlo) were done separately and the anal- 
ysis of each was done by the two authors independently. 
Putting these two unbiased results together, we propose 
a final value of 162.5 ± 2, where the estimated error is 



[1] J. W. Essam, Rep. Prog. Phys. 43, 833 (1980). 

[2] D. StauflFer and A. Aharony, Introduction to Percolation 

Theory (Taylor & Francis, London, 1992), 2 ed. 
[3] V. Privman, P. C. Hohenberg, and A. Aharony, in Phase 

Transitions and Critical Phenomena, edited by C. Domb 

and J. L. Lebowitz, vol. 14, pp. 1-134 (Academic, New 

York, 1991). 

[4] M. Casells, P. Grinza, and A. Rago, J. Stat. Mech.: Th. 

Exp., P10009 (2004). 
[5] G. Delfino, J. Phys. A 37, R45 (2004). 
[6] G. Delfino and P. Grinza, Nucl. Phys. B 682. 521 (2004). 
[7] I. G. Enting and A. J. Guttmann, Physica A 321, 20 

(2003). 

[8] K. A. Seaton, J. Stat. Phys. 107, 1255 (2002). 
[9] L. N. Shchur and O. A. Vasilyev, Phys. Rev. E 65, 016107 
(2001). 

[10] S. B. Lee, Phys. Rev. E 53, 3319 (1996). 

[11] D. Fioravanti, G. Mussardo, and P. Simon, Phys. Rev. E 

63, 016103 (2001). 
[12] G. Delfino and J. L. Cardy, Nucl. Phys. B 519, 551 

(1998). 

[13] M. Corsten, N. Jan, and R. Jerrard, Physica A 156, 781 
(1989). 

[14] H. Kesten, Commun. Math. Phys. 74, 41 (1980). 
[15] M. E. J. Newman and R. M. Ziflt, Phys. Rev. Lett. 85, 
4104 (2000). 

[16] A. J. Guttmann, in Phase Transitions and Critical Phe- 
nomena, edited by C. Domb and J. L. Lebowitz, vol. 13, 
pp. 1-234 (Academic, New York, 1989). 

[17] D. Daboul, A. Aharony, and D. Staufler, J. Phys. A 33, 



5 



1113 (2000). [19] J. L. Cardy and R. M. Ziff, J. Stat. Phys. 110, 1 (2003). 

[18] G. T. Delfino, G. T. Barkema, and J. L. Cardy, Nucl. 
Phys. B 565, 521 (2000). 



